Threshold Effects in Multi-channel Coupling and Spectroscopic Factors in Exotic 

Nuclei 



N. Michelttll, W. Nazarewiczt**, M. Ploszajczak§ 
^ Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996, USA 
* Physics Division, Oak Ridge National Laboratory, P. O. Box 2008, Oak Ridge, TN 37831, USA 
^ Joint Institute for Heavy Ion Research, Oak Ridge, TN 37831, USA 
* Institute of Theoretical Physics, Warsaw University, ul. Hoza 69, 00-681 Warsaw, Poland 
^ Grand Accelerateur National d'lons Lourds (GANIL), 
CEA/DSM - CNRS/IN2P3, BP 55027, F-14076 Caen Cedex, France 
(Dated: February 9, 2008) 

In the threshold region, the cross section and the associated overlap integral obey the Wigner 
threshold law that results in the Wigner-cusp phenomenon. Due to flux conservation, a cusp anomaly 
in one channel manifests itself in other open channels, even if their respective thresholds appear at 
a different energy. The shape of a threshold cusp depends on the orbital angular momentum of 
a scattered particle; hence, studies of Wigner anomalies in weakly bound nuclei with several low- 
lying thresholds can provide valuable spectroscopic information. In this work, we investigate the 
threshold behavior of spectroscopic factors in neutron-rich drip-line nuclei using the Gamow Shell 
Model, which takes into account many-body correlations and the continuum effects. The presence 
of threshold anomalies is demonstrated and the implications for spectroscopic factors are discussed. 
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In 1948, based on general principles, Wigner pre- 
dicted [l| a characteristic behavior (a cusp) of scattering 
and reaction cross sections in the vicinity of a reaction 
threshold. This particular behavior (often referred to 
as the Wigner's threshold law or the Wigner-cusp phe- 
nomenon}_ was given a quantitative explanation a decade 
S B particular, it has been noted 

that, due to the unitarity of the scatter- 
ing matrix and the resulting flux conservation, the pres- 
ence of a threshold anomaly in an opening reaction chan- 
nel can trigger an appearance of Wigner cusps in other 
open channels with lower thresholds. As the shape of the 
cusp strongly depends on orbital angular momentum (is 
strongest in s and p waves), it was early realized that 
the presence of cusp anomaly could provide structural 
information about reaction products [3, 0j ■ 

The Wigner-cusp phenomenon has been studied exper- 
imentally and theoretically in various areas of physics: 
pion- nucleus scattering [§, [l3|, electron-molecule scat- 
tering [ll[, electron-atom scattering [l^], and ultra- 
cold atom-diatom scattering [3]. In low-energy nuclear 
physics, threshold effects have been investigated in, e.g., 
charge-exchange reactions fTJ| , neutron elastic scatte ring 
p^ . and deuteron stripping [16]. Abramovich et al. [17| 
reviewed threshold phenomena in nuclear reactions, in- 
cluding those with the light neutron-rich systems such 
as ^He, i°Be, and i°Li (see also Ref. [ll|) that are of 
particular interest in the context of this work. The in- 
fluence of threshold effects on the cross section and the 
strength function, hence the spectroscopic factor (SF) of 
a threshold state was pointed out in Refs. [§. [Toj. 

The purpose of this study is many-fold. Firstly, we in- 
vestigate whether the Wigner-cusp phenomenon appears 
naturally in a microscopic many-body approach rooted 
in an effective inter-nucleon interaction. The second goal 



is to investigate the influence of threshold effects and the 
coupling to the non-resonant continuum on overlap inte- 
grals, or spectroscopic factors. We demonstrate that the 
energy dependence of SFs caused by an opening of a re- 
action channel can be described at the shell-model level 
only if nucleon-nucleon correlations involving scattering 
states are treated properly. Finally, we emphasize the 
importance of experimental studies of cusp phenomena 
in weakly bound, neutron-rich nuclei in which low-lying 
one-neutron and two-neutron thresholds appear. 

The traditional shell model (SM) of nuclear structure 
views the nucleus as a closed quantum system (CQS) in 
which nucleons occupy bound orbits. While such an as- 
sumption may be somehow justified for well-bound nuclei 
having high particle-emission thresholds, it can no longer 
be applied to weakly bound or unbound systems in the 
vicinity of drip lines, where the coupling to the particle 
continuum (both resonances and the non-resonant scat- 
tering states) becomes important. This coupling can be 
considered in the open quantum system (OQS) extension 
of the SM, the so-called continuum SM (see Ref. [20*] for 
a recent review). In this work, we apply the complex- 
energy implementation of the continuum SM, the so- 
called Gamow Shell Model (GSM) [HI [13 in the ver- 
sion of Refs. [21!, [1^. GSM is a multi-configurational SM 
with a single-particle (s.p.) basis given by the Berggren 
ensemble [24,] which consists of Gamow (bound and res- 
onance) states and the non-resonant continuum. For a 
given Hamiltonian, the number of particles occupying 
states of the non-resonant continuum is a result of GSM 
variational calculations. The resonant states of the GSM 
are the generalized eigenstates of the time-independent 
Schrodinger equation which are regular at the origin and 
satisfy purely outgoing boundary conditions. The GSM 
can thus be viewed as a quasi-stationary many-body OQS 
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formalism. (An alternative microscopic approach, suc- 
cessfully applied to the structure of weakly bound or un- 
bound nuclei, is the Complex Scaling Method that is ca- 
pable of treating different kinds of reaction channels and 
continuum states starting from different thresholds [1^.) 
Since the Wigner-cusp phenomenon is most pronounced 
for low-^ waves and for neutrons (no Coulomb barrier), 
as an illustrative example we choose the case of the 
one-neutron (In) channel in the model two- and three- 
neutron systems outside the inert core: ^He and ''He. 
Our aim is not to fit the actual experimental data, but 
rather to accomplish the physics goals as stated above. 

Gamow Shell Model Framework- The starting point 
of GSM is the Berggren one-body completeness relation 
[23 | allowing the expansion of both bound and unbound 
states. The Berggren ensemble consists of resonant and 
scattering states generated by a finite-depth potential 
V{r). Resonant states are solutions of the Schrodingcr 
equation with purely outgoing asymptotics. Their en- 
ergies and wave functions are in general complex. The 
resonant states of the GSM have either bound or de- 
caying character; they form the so-called pole subspace. 
Scattering states entering the Berggren ensemble are also 
defined in the complex energy/momentum plane. For a 
given partial wave {£,j), the scattering states are dis- 

£ ■ 

tributed along a contour in the complex momentum 
plane. The set of all bound and decaying states |u„) en- 

closed between L_j^ and the real fc-axis, and scattering 

states \uk) in is complete j2^ : 

\ub){^\ = 1, (1) 

where \B) is either a discrete resonant state or a scatter- 
ing continuum state (L^ part). The tilde symbol above 
bra vectors in Eq. ([T]) signifies that the complex conju- 
gation arising in the dual space affects only the angular 
part and leaves the radial part unchanged [24| . The con- 
tinuous part of the completeness relation ([T]) has to be 
discretized in numerical applications. For that purpose, 
scattering states are discretized and renormalized in or- 
der to obtain a discrete completeness relation Q : 

JV 

; Wb) = Wke), (2) 

where {fcg, cjg} is the set of discretized momenta and as- 
sociated weights provided by a Gauss-Legendre quadra- 
ture. The many-body GSM basis corresponds to Slater 
determinants (SD) spanned by one-body Berggren states: 
\SDi) = \ui-^ ■ ■ -Ui^) where \SDi) is the i-th SD in the 
A-body basis and Ui is the j-th one-body state occupied 
in 15*1?^). The many-body completeness relation is built 
from Eq. ^ by forming all possible SDs generated by 
the Gamow one-body states: 

5]|5A)(5A|^1. (3) 



The completeness in ^ is not exact as the one-body 
completeness relation ^ is approximate due to the dis- 
cretization. In the basis Q, the GSM Hamiltonian H 
becomes a complex symmetric matrix. Moreover, many- 
body bound and resonant states are embedded in the 
background of non-resonant scattering eigenstates, so 
that one needs a criterion to isolate them. The overlap 
method [2l| has proven to be very efficient to solve this 
problem. For this, one diagonalizes first H in the pole 
subspace to generate a zeroth-order vector Ivto)- In the 
second step, |5'o) is used as a pivot to generate a Lanc- 
zos subspace of the full GSM space. Its diagonalization 
provides eigenvectors of H in the total GSM space, and 
the requested bound or decaying eigenstate of H is the 
one which maximizes the overlap | (^'o 1 5*) | ■ 

The definition of observables in GSM follows directly 
from the mathematical setting of quantum mechanics in 
the rigged Hilbert space rather than the usual Hilbert 
space [26;,^]. Modified definition of the dual space, em- 
bodied by the tilde symbol above bra states, implies that 
observables in many-body resonances become complex. 
In this case, the real part of a matrix element corresponds 
to the expectation value, and the imaginary part can be 
interpreted as the uncertainty in the determination of 
this expectation value due to the possibilityof decay of 
the state during the measuring process [24|.|28|. 

GSM Hamiltonian- In this Letter, the s.p. basis ([!]) 
is generated by a Woods-Saxon (WS) potential with the 
radius Ro—2 fm, the diffuseness d=0.65 fm, the spin-orbit 
strength Vso=7.5 MeV, and the depth of the central po- 
tential Vo=47 MeV (the "^He" parameter set). This po- 
tential reproduces experimental energies and widths of 
the s.p. resonances 3/2^ and l/2;f in ^He. The GSM 
Hamiltonian is a sum of the one-body WS potential, 
representing the effect of an inert ^He core, and of the 
two-body interaction among valence particles, given by a 
finite-range surface Gaussian interaction (SGI) [11] with 
the range /i = 1 fm and the coupling constants depend- 
ing on the total angular momentum J of the neutron 
pair: V^°^ = -403 MeV fm^ and V^^^ = -392 MeV 
fm'^. These constants are fitted to reproduce the exper- 
imental ground state (g.s.) binding energies of ^He and 
^He with the "^He" WS potential. The valence space 
for neutrons consists of the P3/2 and pi/2 partial waves. 
The wave functions include a OP3/2 resonant state 
and P3/2 non-resonant scattering states along a complex 
contour enclosing the OP3/2 resonance in the complex k- 
plane. For a pi/2 part, we take non-resonant scattering 
states along the real- A: axis (the broad Opi/2 resonant 
state plays a negligible role in the g.s. wave function of 
^He). For both contours, the maximal momentum value 
is kniax =3.27 fm~^. The contours have been discretized 
with up to 60 points and the attained precision on ener- 
gies and widths is better than 0.1 keV. In the subsequent 
analysis, parameters of the GSM Hamiltonian are varied 
in order to change positions of In thresholds in various 
isotopes. 

In order to illuminate the continuum coupling effects in 
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GSM, we have introduced a simphfied harmonic oscillator 
SM (HO-SM) scheme. Here, the radial wave functions 
are those of the spherical harmonic oscillator with the 
frequency huj — 41A'~^^'^ MeV. The one-body part of 
the HO-SM Hamiltonian is given by the real energies of 
the one-body part of the GSM Hamiltonian. The HO- 
SM scheme is supposed to illustrate the "standard" CQS 
SM calculations in which only bound valence shells are 
considered in the s.p. basis. 

Spectroscopic Factors in GSM- SFs are useful indica- 
tors of the configuration mixing in the many-body wave 
function. Extensive attempts have been made to deduce 
SFs using direct reactions, such as single- nucleon trans- 
fer, nucleon knockout, and elastic break-up reactions, us- 
ing hadronic and leptonic probes. These analyses often 
reveal model- and probe-dependence [l^, [13, |3l[ raising 
concerns about the possibility of their precise experimen- 
tal determination. (For an extensive study of spectro- 
scopic factors in exotic nuclei from nucleon-knockout re- 
actions, see Ref. 32].) Studies of (e, e'p) reactions in 
closed-shell nuclei [33] demonstrated that the SFs are re- 
duced by ~ 35% with respect to the standard SM pre- 
dictions, mainly due to the coupling to high-momentum 
states reached by the short-range and tensor components 
of the nucleon-nucleon interaction [s^ . In this work, we 
point out that additional difficulties in extracting and 
interpreting SFs from the measured cross sections using 
the standard SM lie in the neglect of particle continuum, 
channel coupling, and strength fragmentation. 

Single-nucleon overlap integrals and the associated 
spectroscopic factors (SFs) are basic ingredients of the 
theory of direct reactions (single-nucleon transfer, nu- 
cleon knockout, elastic break-up) (ssl . [36j . Experimen- 
tally, SFs can be deduced from measured cross sec- 
tions; they are useful measures of the configuration mix- 
ing in the many-body wave function. The associated 
reaction-theoretical analysis often reveals model- and 
probe-dependence [l^, [Sfl, HH raising concerns about 
the accuracy of experimental determination of SFs. In 
our study we discuss the uncertainty in determining SFs 
due to the two assumptions commonly used in the stan- 
dard SM studies, namely (i) that a nucleon is transferred 
to/from a specific s.p. orbit (corresponding to an ob- 
served s.p. state), and (ii) that the transfer to/from the 
continuum of non-resonant scattering states can be dis- 
regarded. In this work, we define SFs in a usual way, 
through the radial overlap functions uej (r) [H, [13, [11] : 



channel function in the complete Berggren basis, one ob- 
tains: 



-A-1 



Ja 



(4) 



where l^*;^-*) 



and I'^A^^i) are wave functions of nuclei 



A and A ^ 1, respectively, and is the angular-spin 
part of the channel function. The angular-spin degrees of 
freedom are integrated out in Eq. ^ so that ugj depends 
only on the relative radial coordinate of the transferred 
particle r = ]7^. The spectroscopic factor, denoted by 
S'^, is defined as usual through the norm of the overlap 
function [1^[13, [11] . Using a decomposition of the 
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where a^j{B) is a creation operator associated with a s.p. 
Berggren state \B). Since Eqs. (|5l6p involve summation 
over all discrete Gamow states and integration over all 
scattering states along the contour L^, the final result 
is independent of the s.p. basis assumed. This is in con- 
trast to standard SF experimental extraction and SM 
calculations where model-dependence enters through the 
specific choice of a s.p. state a^^-, with Eq. © reducing 

to the sole matrix element {'^^/ ] ]a+^^. W^''/^^)"^ . This can 
lead to sizeable errors if the states of ^ — 1 and/or A lie 
close to a channel threshold. 

The Threshold Cusp- The GSM results for *^He 
g.s. SF in the channel: J^He(g.s.) ® ^3/2]*^^ are shown 
in Fig. [IJi as a function of one-neutron separation energy 
Sin in ^He. To this end, we have fixed the depth of the 
WS potential so that 0^3/2 and Opi/2 are both bound 
with respective energies -5 and -0.255 MeV, and we have 
varied the SGI monopole coupling constant Vy~'^^ so 
that the In separation energy of ^He goes through zero. 
At S'in=0 (In-emission threshold), the calculated SF 
exhibits behavior consistent with the Wigner threshold 
law: the quickly varying component of SF behaves as 
{—SinY^^^'^ (^=1) below the In threshold and follows the 
{SinY'^^^'^ rule above the threshold. It is worth noting 
that the parameters of the GSM Hamiltonian and the as- 
sociated S-matrix poles do not show any discontinuities 
around S'i„=0. This result constitutes an excellent test of 
the GSM formalism: the Wigner limit is reached precisely 
at the In threshold obtained from many-body calcula- 
tions. It is interesting to see that the HO-SM results are 
depressed by about 25% as compared to GSM. Indeed, in 
the HO-SM s.p. basis, the configurations ]0p3/2'8)0p3/2]°^ 

and ]0pi/2'8)0pi/2]°^ are strongly mixed by the SGI resid- 
ual interaction; hence, the value of p^/2 SF in the g.s. of 
^He is significantly reduced. Moreover, the s.p. basis in 
HO-SM calculations comprises HO states whose radial 
form factors are independent of the depth of the WS po- 
tential. Hence, no threshold effect can be seen in HO-SM 
SFs. 

Threshold Effects Due to Channel Coupling- As dis- 
cussed above, when a new channel opens at energy Et^ 
there appears a flux redistribution in other open channels 
with lower thresholds. This flux redistribution may be 
affected by the presence of the Wigner cusp in the new 
channel; hence, the reaction cross-sections in all open 
channels may exhibit the threshold anomaly at Et- To 
illustrate the phenomenon of channel coupling. Fig. 
shows again the SF for the J^He(g.s.) ^^3/2]''^ channel 
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FIG. 1: The real part of the overlap integral as a func- 
tion of one-neutron separation energy S\n as indicated. 

Top and middle: (''He(g.s.)|[^He(g.s.) ® P3/2]''^>; bottom; 

(^He(g.s.)|[®He(g.s.) (g) p3/2]°^>. Solid line: GSM results; 
dashed line: the SM-like approximation (HO-SM) where 
the SGI matrix elements are calculated in the HO basis 
{0p3/2, Opi/2}. The In thresholds in ^He (top and bottom) 
and ^He are marked by arrows. The results displayed in 
the middle panel are plotted as a function of the energy of 
the 0p3/2 resonant state, which, in our model, is negative of 
Si„[^He]. For more details, see the discussion in the text. 



but this time as a function of the OP3/2 s.p. state energy 
(or negative S\n of ^He). This calculation was carried 
out by varying the depth of the WS potential so that the 
P3/2 pole of the S'-niatrix (which is also the g.s. of ^He 
in our model space), would change its character from a 
bound state to an unbound decaying Gamow state. At 
'S'in[^He]=0, the In and 2n thresholds in ^He become de- 
generate, i.e., the two channels of ^He-|-n, namely ^He 
and ^He-|-2n couple. As seen in Fig. [1] the SF in the 
pHe(g.s.) (8'P3/2]° channel exhibits the Wigner cusp at 
S'i„[^He]=0, i.e., at a 2n threshold. At this point, the 
first derivative dS'^ / de^p^^^ of the SF becomes infinite, 
consistent with the Wigner law for £=1. The SF in HO- 
SM varies very little in the whole studied range of 0^3/2 
energies. (Let us mention that the competition between 
bound and unbound states in ^He, and the ^He-fn and 
^He-f2n reaction channels has been studied within the 



coupled cluster method |25|.) 

In the case shown in Fig. [T|d the GSM Hamiltonian 
changes at S'i„[^He]=0 as at this point the OP3/2 pole 
becomes unbound. It is, therefore, instructive to inves- 
tigate the situation at which the Hamiltonian behaves 
smoothly around the threshold. Figure (TJ; illustrates the 
case of ''He g.s. in the channel [^He(g.s.) ® Pz/2f'^'^ ■ 
In order to avoid secondary open-channel mixing effects, 
we have adjusted the WS potential depth so that the 
^He g.s. is bound by 5 MeV and cannot play any role 
in the anomalous energy dependence of the studied SF. 
The Opi/2 state is weakly bound with an energy of —0.255 
MeV. In order to control the binding energy of ^He, the 
^(-'-n) coupling constant is varied so that ^He g.s. is 
always bound while Sin of ^He goes through zero. The 
GSM space is the same as in the previous cases, except 
that the Qpi/2 state is now included in the GSM basis as it 
is bound. A cusp in the calculated SF for ^He is clearly 
seen at S'i„[^He]=0. The threshold anomaly shown in 
Fig. [TJ; can only result from the cross-channel couplings. 
Again, the SF in HO-SM smoothly varies in the whole 
energy region considered. 

Conclusions- The anomalies in the spectroscopic fac- 
tors when the total energy of the system varies through 
the threshold of an opening channel are discussed within 
the many-body OQS formalism of the GSM. The main 
conclusions of this work can be summarized as fol- 
lows: (i) Many-body OQS calculations correctly predict 
the Wigner-cusp and channel-coupling threshold effects. 
This constitutes a very strong theoretical check for the 
GSM approach; (ii) The spectroscopic factors defined in 
the OQS framework through the norm of the overlap in- 
tegral, exhibit strong variations around particle thresh- 
olds. Such variations cannot be described in a standard 
CQS SM framework that applies a 'one-isolated-state' 
ansatz and ignores the coupling to the decay and scatter- 
ing channels. In our model calculations, the contribution 
to SF from a non-resonant continuum can be as large as 
25%; (iii) Any theoretical model aiming at a meaning- 
ful description of SFs of \ow-£ states in weakly bound 
nuclei must meet certain minimal conditions. Namely, 
it should be able to account for (a) many-body config- 
uration mixing and the resulting spreading of the spec- 
troscopic strength due to inter-nucleon correlations, (b) 
coupling to the particle continuum that affects radial 
properties of wave functions in the neighborhood of re- 
action thresholds, and (c) coupling between various re- 
action channels. 

Considering points (ii) and (iii) above, one should be 
careful when extracting spectroscopic information from 
transfer experiments on drip-line nuclei. The results pre- 
sented in this work suggest that, similar to other fields, 
experimental studies of various aspects of threshold ef- 
fects could provide valuable spectroscopic information 
about the s.p. structure of weakly bound nuclei. Here, of 
particular interest are systems with near-lying In and 2n 
thresholds such as ®^®He or non-Borromean two-neutron 
halos [39] , in which cusps in SFs are expected to be par- 
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ticularly strong in \ow-£ {£ = 0, 1) neutron channels. Sim- 
ilar features were found in the analysis of the continuum- 
coupling correction to the CQS eigenvalues near the reac- 
tion threshold [i^]. Finally, let us note that a threshold 
anomaly is also expected in proton-rich nuclei. While the 
effect is weaker than in the neutron-rich systems, one still 
expects anomalous differences in SFs and other spectro- 
scopic quantities in mirror nuclei. Work along these lines 



is in progress. 
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